Inclinations and black hole masses of Seyfert 1 galaxies 



Xue-Bing Wu 1 and J. L. Han 2 

1. Department of Astronomy, Peking University, Beijing 100871, China; 

wuxb@bac.pku. edu. cn 

2. National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, 

China; hjl@bao.ac.cn 

ABSTRACT 

A tight correlation of black hole mass (Mbh) and central velocity dispersion 
(a) has been found recently for both active and quiescent galaxies. By applying 
this correlation, we develop a simple method to derive the inclination angles for a 
sample of 11 Seyfert 1 galaxies that have both measured central velocity disper- 
sions and black hole masses estimated by reverberation mapping. These angles, 
with a mean value of 36° that agrees well with the result obtained by fitting 
the iron Ka lines of Seyfert Is observed with ASCA, provide further support to 
the orientation-dependent unification scheme of AGN. A positive correlation of 
the inclinations with observed FWHMs of B./3 line and a possible anti-correlation 
with the nuclear radio-loudness have been found. We conclude that more ac- 
curate knowledge on inclinations and broad line region dynamics is needed to 
improve the black hole mass determination of AGN with the reverberation map- 
ping technique. 

Subject headings: black hole physics - galaxies: active - galaxies: nuclei - galax- 
ies: Seyfert. 



1. Introduction 

Supermassive black hole is believed to be an essential part of active galaxy and quasar 
(Lynden-Bell 1969; Rees 1984). Recently a lot of evidence has also been found for the 
existence of supermassive black holes in nearby quiescent galaxies (Kormendy & Richstone 
1995; Magorrian et al. 1998). A significantly tight correlation of black hole mass (Mbh) 
with the bulge velocity dispersion (a) was found for nearby galaxies (Gebhardt et al. 2000a; 
Ferrarese & Merritt 2000). More recent studies indicated that some Seyfert galaxies with 
Mbh measured by reverberation mapping method follow the same Mbh-c relation as for 
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normal galaxies (Gebhardt et al. 2000b; Ferrarese et al. 2001), which implies that this 
relation may be universal for galaxies. This relation strongly suggests a tight connection 
between the formation and evolution of the supermassive black hole and the galactic bulge, 
though the nature of this connection is still unclear. 

The black hole masses of some active galactic nuclei (AGN) have been recently estimated 
by the reverberation mapping technique (Wandel, Peterson & Malkan 1999; Ho 1999; Kaspi 
et al. 2000). With this technique, the broad line region (BLR) size can be measured using 
the lag between the variability of continuum and emission line fluxes. The black hole mass 
can be then estimated from the BLR size and the characteristic velocity (determined by the 
full width at half-maximum (FWHM) of emission line). However, it may not be so strait- 
forward to estimate the characteristic velocity of BLR according to the observed FWHMs 
of emission lines (Fromerth & Melia 2000). Many effects, especially the inclination and 
BLR geometry, can lead to larger uncertainties to the estimation of black hole mass using 
reverberation mapping (Krolik 2001). 

Inclination is an important ingredient in the orientation dependent unified scheme of 
AGN (Urry & Padovani 1995), but it is not easy to be derived directly from observations. 
By fitting the observed iron Ka line profile with the accretion disk model, Nandra et 
al. (1997) estimated the inclinations, with a mean value of 30°, for 18 Seyfert 1 galaxies. 
Although with large uncertainties, their result suggested that the inclinations of Seyfert Is are 
small, consistent with Seyfert 1/2 unification scheme (Antonucci & Miller 1985). Moreover, 
numerous observations have shown the evidence for the lack of edge-on Seyfert 1 galaxies 
(Keel 1980) and suggested that Seyfert 1 galaxies have the dusty torus with half opening 
angles of 40° to 60° (Osterbrock & Martel 1993; Ho, Fillippenko & Sargent 1997; Schmitt 
et al. 2001). Therefore, the inclinations of Seyfert 1 galaxies are generally expected to be 
small, though further evidence is obviously needed. 

If the Mbh-0" relation is valid for Seyfert galaxies, this provides an independent way 
to determine Mbh from the reverberation mapping. Using Mbh derived from the measured 
stellar velocity dispersion and assuming a simple BLR geometry, we develop a method to de- 
rive the inclinations for 11 Seyfert 1 galaxies with both measured central velocity dispersions 
and estimated black hole masses by reverberation mapping. Our results are supported by 
several correlation studies and are consistent with previous knowledge of inclination effects 
of AGN. 
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2. Black hole mass determinations of AGN 

According to the reverberation mapping technique, the black hole mass can be estimated 
by a virial form, 

V 2 R 

M rev = (1) 

provided that the BLR region of AGN is gravitational bounded and has a Keplerian char- 
acteristic velocity (V). The BLR size (R) can be derived by the lag between the variability 
of continuum and broad emission line. Assuming AGN having random inclinations, Wandel 
et al. (1999) and Kaspi et al. (2000) related the BLR characteristic velocity to the FWHM 
of H/3 emission line by V = (\/3/2)Vfwhm- However, as pointed recently by McLure & 
Dunlop (2001), the assumption of random inclinations seems unrealistic for quasars. Many 
observational studies also indicated that Seyfert 1 galaxies are probably not viewed at all 
inclinations (Osterbrock & Martel 1993; Ho et al. 1997; Nandra et al. 1997; Schmitt et 
al. 2001). 

On the other hand, the BLR dynamics of AGN has not been well understood yet. The 
simplest case may be described as circular orbits confined to the disk plane, but the real 
case is probably much more complicated. The same as that assumed for quasars, the BLR 
velocity of Seyfert Is may be better represented by a combination of a random isotropic 
component, with characteristic velocity V r , and a component in the disk plane, with char- 
acteristic velocity V p (Wills & Browne 1986). Therefore, the observed FWHM of emission 
line will be given by 

^FWHM = 2(V; 2 + V; 2 sin^) 1 / 2 , (2) 

where % is the inclination angle of the disk normal relative to the line of sight. Assuming 
that V p is significantly larger than V r and that % lies randomly between and 46°, McLure 
& Dunlop (2001) has reproduced the distribution of observed H/3 FWHMs for a sample of 
AGN using the above formula. If we relate the observed emission line FWHM and BLR 
characteristic velocity with eq. (2), the black hole mass of AGN can be obtained by 

4(sm i + A z ) G 

where A = V r /V p . In most recent studies, the black hole mass (M rev ) determined by re- 
verberation mapping is based on the assumption of random inclinations and is expressed 
as M rcv = l ^wjpg (Wandel et al. 1999; Kaspi et al. 2000). If the inclination effect is 
considered, using eq. (3) we can derive the inclination angle by: 



arcsin( A /^--A2). (4) 
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If M B h of AGN can be derived independently by other methods, we can use eq. (4) to 
estimate the inclinations of AGN. Fortunately a tight relation between black hole mass and 
central velocity dispersion seems to exist for AGN, therefore Mbh can be directly derived 
from the measured central velocity dispersion. According to Gebhardt et al. (2000a), the 
Mbh-0" relation is: 

M BH = 1-2 x 1O 8 M (ct/2OO km s^ 1 ) 3 ' 75 . (5) 

Using a slightly steeper slope found by Ferrarese & Merritt (2000) has no significant effect 
on our results. 



3. Inclination angles of Seyfert 1 galxies 

We collected 11 Seyfert 1 galaxies with both measured central velocity dispersion (a) 
and estimated black hole masses (M rev ) by reverberation mapping (Wandel et al. 1999; Ho 
1999; Kaspi et al. 2000). Their a and M rcv values, together with the observed FWHMs of 
H/9 line (Wandel et al. 1999; Ho 1999) and the nuclear radio-loudness (Ho & Peng 2001), 
are summarized in Table 1. Seven sources have a values from high quality measurements 
of Ferrarese et al. (2001) and Di Nella et al. (1995). The a values of other four sources 
were adopted from Nelson & Whittle (1995). The uncertainties of M rev and V F whm for 
NGC 3516 and NGC 4593 are unavailable in literature and were assumed to be 30% and 
5%, respectively. 

The inclinations of these 11 Seyfert Is can be estimated using eq. (4) and eq. (5). 
For simplicity we further assume A <C M rev /3M B H in eq. (4), which is equivalent to the 
approximation that the BLR characteristic velocity is dominated by the component in the 
disk plane (McLure & Dunlop 2001). Under this assumption, we derived the inclinations 
for 11 Seyfert Is (see Table 1 and Figure 1). The errors of i were estimated from the 
uncertainties of both o and M rev . The inclination angles we derived are in a range from 
20° to 60°, with a mean value of 36.2°. This agrees with the result obtained by fitting the 
observed iron Ka line in X-ray band with accretion disk model (Nandra et al. 1997), and 
is consistent with the expectation of unified scheme of AGN (Antonucci & Miller 1985). In 
Figure 2, we can clearly observe an apparent positive correlation between inclinations and 
observed H/3 FWHMs. A minimum \ 2 fit considering the errors of both parameters gives, 
sin(i) = (0.23 ± 0.09) + (0.086 ± 0.024) (V FWH M/1000km s" 1 ), with X 2 and probability of 
7.13 and 62%, respectively. The zero slope is less favored because it produces a larger x 2 
(22.78) with a probability of 1% only. A simple Spearman test gives the correlation coefficient 
R = 0.78, with have a probability of P = 4.47x 10~ 3 to occur by chance. We have also applied 
the bootstrap method to investigate the uncertainty of correlation coefficient by considering 
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the uncertainties of both parameters, and obtained (R) = 0.58 ± 0.18, which indicates 
a moderately significant correlation. In Figure 3, we plotted the inclinations against the 
nuclear radio loudness, defined by the ratio between 5 GHz nuclear radio luminosity and B- 
band nuclear optical luminosity (Ho & Peng 2001). Because there are serious contaminations 
to the luminosity of Seyfert nucleus from the host galaxy, using the nuclear radio-loudness 
can maximumly diminish such contaminations and better describe the nature of Seyfert 
nuclei. Although only seven sources with available data for nuclear radio-loudness, Figure 3 
still displays the trend that Seyfert Is with larger nuclear radio-loudness may have smaller 
inclinations. These results seem to be consistent with our knowledge on the inclination 
effects of quasars and AGN (Wills & Browne 1986; Wills & Brotherton 1995). 

If the characteristic velocity of BLR is dominated by the component in the disk plane, 
the intrinsic FWHM of H/3 emission line can be approximated by the observed FWHM 
divided by sin(i). In Figure 4 we show the distributions of intrinsic FWHMs of H/3 line and 
M BH . Two narrow line Seyfert Is (NLSls), NGC 4051 and Mrk 110, locate in the lower-left 
corner of Figure 4, though their inclination angles, estimated to be 19.6° and 37.4°, are 
not significantly smaller than those of other broad line Seyfert Is. This indicates that the 
narrowness of emission lines of NLS1 may not be simply due to the inclination effects, but is 
probably more related to their smaller black hole masses. These NLSls may have accretion 
rates close to Eddington limit (Boiler, Brandt & Fink 1996). For most broad line Seyfert Is, 
the Eddington ratios (L/ L E dd) are in the range of 0.1 to 0.01 (Wandel et al. 1999), which 
are significantly lower than those of NLSls. 



4. Discussion 

Assuming a Kelperian rotation of BLR and Seyfert Is following the same Mbh-o" relation 
as normal galaxies, we develop a method to derive the inclinations of 11 well studied Seyfert 
1 galaxies. The values of these inclinations derived by us seem to be supported by a positive 
correlation with the observed H/5 FWHMs and a possible anti-correlation with the nuclear 
radio loudness. Our results agree well with previous knowledge on inclinations of Seyfert Is 
and are consistent with the expectation of unification scheme of AGN. 

It is worth noting that the BLR dynamics may not be as simple as we assumed, though 
the Keplerian rotation of BLR have been confirmed by some recent studies (Peterson & 
Wandel 1999, 2000). The inclinations we derived for 11 Seyfert Is may be regarded as up- 
per limits because we assumed the characteristic velocity in the disk plane is the dominant 
component in BLR. Considering of the isotropic component will reduce the values of incli- 
nations (see eq. (4)). VLBI study on the radio- loud Seyfert 1 galaxy 3C 120 derived a lower 
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inclination angle (about 10°) based on the synchrotron self-Compton jet model (Ghisellini 
et al. 1993), which is smaller than but still marginally consistent with the value obtained by 
us. 

We noticed that ratio between the black hole mass determined by reverberation mapping 
(M rev ) and the 'true' black hole mass (M B h), can be approximated by 3sin 2 (i) (see eq. (4)). 
With the mean value of i (36.2°) derived by us, we obtained M rev /MBH — 1-05 for Seyfert 
Is. It means that the black hole mass estimated by the standard method of reverberation 
mapping, can still represent the 'true' black hole mass well if the inclination of Seyfert galaxy 
is not substantially different from 36.2°. This may be the reason why Seyfert galaxies also 
follow the Mbh-0" relation even if the values of Mbh measured by reverberation mapping 
were adopted (Gebhardt et al. 2000b; Ferrarese et al. 2001). However, the M rev /M BH 
value changes from 0.35 to 2.25 if the inclination increases from 20° to 60°. Therefore, more 
accurate information about inclinations and the BLR dynamics of AGN will be undoubtedly 
helpful to the improvement of central black hole mass estimations with the reverberation 
mapping technique. 
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Science Foundation, and the National Key Basic Research Science Foundation (NKBRSF 
G19990752) in China. 
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Table 1: Data of Seyfert 1 galaxies 

Name a Ref* M rev Vfwhm(^^s) log R nuc Inclination 



(km/s) (1O 7 M ) (km/s) (°) 



3C 120 


162±20 


1 


2-311:5 


2210±120 




22.01?:? 


Mrk 79 


130±9 


2 


c O+2.0 

O . Li 2 g 


6280±850 




58.5l27:g 


Mrk 110 


86±5 


2 


0.561°-? 


1670±120 




37.41^ 


MrK 590 


169±28 


1 


l-78±g : lf 


2170±120 


1.62 


17-81^ 


Mrk 817 


142±6 


2 


4.41}-? 


4010±180 


1.21 


4i.6i?:l 


NGC 3227 


144±22 


1 


<->• z> _3 g 


5530±490 


1.12 


q 7 r+17.3 
° ' - d -25.4 


NGC 3516 


124±5 


3 


9 q+0.69 
z -°-0.69 


4760±240 


0.78 


38.3±7.6 


NGC 4051 


80±4 


2 


o n+ 013 


1230±60 


0.87 


19.6™ 


NGC 4151 


93±5 


2 


i cq+1.06 


5230±920 


0.49 


6o.oii8:8 


NGC 4593 


124±29 


1 


0.811°- 


3720±180 




21.6±10.5 


NGC 5548 


183±10 


2 


12.3lfJ 


5500±400 


1.24 


^q 7 +7.6 

'-6.9 



* References: 1, Nelson & Whittle (1995); 2, Ferrarese et al. (2001); 3, Di Nella et al. 
(1995). 
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Fig. 1. — The derived inclination angles of 11 Seyfert 1 galaxies against the black hole masses 
determined by the M B1 i-a relation. The dashed line represents the mean value (i) = 36.2°. 
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Fig. 2. — The inclination angles against the observed FWHMs of H/3 line. The dashed line 
shows the minimum x 2 fit with errors of both parameters considered. 
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Fig. 3. — The inclination angles against the nuclear radio-loudness. 
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Fig. 4. — The intrinsic FWHMs of H/3 line against black hole masses determined by the 
Mbh-0" relation. 



